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Abstract ^ Gilbert tessellation arises by letting linear segments (cracks) in M^ unfold 
in time with constant speed, starting from a homogeneous Poisson point process of germs 
in randomly chosen directions. Whenever a growing edge hits an already existing one, it 
stops growing in this direction. The resulting process tessellates the plane. The purpose of 
i_i the present paper is to establish law of large numbers, variance asymptotics and a central 

limit theorem for geometric functional of such tessellations. The main tool applied is the 
>■ stabilization theory for geometric functional. 
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1 Introduction and main results 

Let A' C M^ be a finite point set. Each x & X is independently marked with a unit length 
random vector cxx making a uniformly distributed angle a^ € [0, tt) with the x-axis, which is 
referred to as the usual marking in the sequel. The collection X = {[x, Oix)}x^x determines 
a crack growth process (tessellation) according to the following rules. Initially, at the time 
t = 0, the growth process consists of the points (seeds) in X . Subsequently, each point 
X E X gives rise to two segments growing linearly at constant unit rate in the directions of 
a^ and — d^ from x. Thus, prior to any collisions, by the time t > the seed has developed 
into the edge with endpoints x — toix and x + td^, consisting of two segments, say the upper 
one [x, X + ta,j\ and the lower one [x, x — ta^]. Whenever a growing segment is blocked by 
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an existing edge, it stops growing in that direction, without affecting the behaviour of the 
second constituent segment though. Since the possible number of collisions is bounded, 
eventually we obtain a tessellation of the plane. The resulting random tessellation process 
is variously called the Gilbert model/tessellation, the crack growth process, the crack 
tessellation, and the random crack network, see e.g. [8, 12] and the references therein. 

Let G{X) denote the tessellation determined by X. We shall write C,~^{x,X), x E X, 
for the total length covered by the upper segment emanating from x in G{X), and likewise 
we let ^~{x,X) stand for the total length of the lower segment from x. Note that we 
use X for marked version of x, according to our general convention of putting bars over 
marked objects. For future use we adopt the convention that if x does not belong to X, 
we extend the definition of ^'^^^(x, X) by adding x to A? and endowing it with a mark 
drawn according to the usual rules. Observe that for some x the values of .^^^^ may be 
infinite. However, in most cases in the sequel X will be a realization of the homogeneous 
Poisson point process V = Vr oi intensity r > in growing windows of the plane. We 
shall use the so-called stabilization property of the functionals C,^ and ^~, as discussed in 
detail below, to show that the construction of G{X) above can be extended to the whole 
plane yielding a well defined process G{V), where, as usual, V stands for a version of V 
marked as described above. This yields well defined and a.s. finite whole-plane functionals 

e+(-,P)andr(-,^). 

The conceptually somewhat similar growth process whereby seeds are the realization 
of a time marked Poisson point process in an expanding window of M^ and which subse- 
quently grow radially in all directions until meeting another such growing seed, has received 
considerable attention [1, 3, 4, 5, 6, 11, 16], where it has been shown that the number of 
seeds satisfies a law of large numbers and central limit theorem as the window size in- 
creases. In this paper we wish to prove analogous limit results for natural functionals 
(total edge length, sum of power- weighted edge lengths, number of cracks with lengths ex- 
ceeding a given threshold etc.) of the crack tessellation process defined by Poisson points 
in expanding windows of M^. We will formulate this theory in terms of random measures 
keeping track not only the cumulative values of the afore-mentioned functionals but also 
their spatial profiles. 

Another interesting class of model bearing conceptual resemblance to Gilbert tessella- 
tions are the so-called lilypond models which have recently attracted considerable attention 
[2, 7, 9, 10] and where the entire (rather than just directional) growth is blocked upon a 
collision of a growing object (a ball, a segment etc.) with another one. 

To proceed, consider a function (j) : [M+UJ-l-oo}]^ — )■ M with at most polynomial growth, 
i.e. for some < q < +oo 

<j){n,r2) = 0{{n+r2y). (1) 

With Qx := [0, vA]^ standing for the square of area A in M^, we consider the empirical 
measure 

^^t■■= E ^{t{x,V),C{x,V))5.^,^^. (2) 



Thus, /i^ is a random (signed) measure on [0, 1]^ for all A > 0. The large A asymptotics of 
these measures is the principal object of study in this paper. Recalling that r stands for 
the intensity of P = P,-, we define 

E(r):=E0(e+(O,P),r(O,P)). (3) 

The first main result of this paper is the following law of large numbers 
Theorem 1 For any continuous function f : [0, 1]^ -^M. we have 

hm - / fdfil = tE{t) / f{x)dx 

in L^, p ^ 1. 

Note that this theorem can be interpreted as stating that E{t) is the asymptotic mass per 
point in /x^, since the expected cardinality oiVnQx is rX. To characterize the second order 
asymptotics of random measures /i^ we consider the pair-correlation functions 

c4x]:=E<p'{e{x,V),r{x,V)), xeR' (4) 

and 

c4x, y] := E0 {^ix, V U {y}).C{x. V U {y})) 
■ {t{y, V u {x}), r(i/, V U {x})) - \E{t)\\ (5) 

In fact, it easily follows by translation invariance that c^\x\ above does not depend on x 
whereas c^fx, y\ only depends owy — x. In terms of these functions we define the asymptotic 
variance per point 



V{t) = c^[0]+r / c^\^,x\d. 

7r2 



x. (6) 

Notice that in a special case when function 0(-, •) is homogeneous of degree k (i.e. for c G M 
we have (f){cri,cr2) = c^(f){ri,r2)) one can simplify (3) and (6). Then the following remark 
is a direct consequence of standard scaling properties of Gilbert's tessellation construction 
and those of homogeneous Poisson point processes, whereby upon multiplying the intensity 
parameter r by some factor p we get all lengths in G{V) re-scaled by factor p~^/^. 

Remark 1 For : [IR+ U {+oo}]^ — ?> M homogeneous of degree k we have 

E{t) = r-'=/2^(l) 
nr)=r-V(l). (7) 

In other words, E{-) and V{-) are homogeneous of degree —k/2 and —k, respectively. 
Our second theorem gives the variance asymptotics for fiy 



Theorem 2 The integral in (6) converges and V{t) > for all r > 0. Moreover, for each 
continuous f : [0, 1]^ — ?> IR 
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tV{t) I f{x)dx. 



Our final result is the central limit theorem 

Theorem 3 For each continuous / : [0, 1]^ — ?■ M the family of random variables 
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converges in law to J\f iO, tV{t) L ^-,2 f'^{x)dx j as A — ?■ oo. Even more, we have 



sup 
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J [0,1] 2 jdl^x 
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for all X > 1, where C is a finite constant. 

Principal examples of functional where the above theory applies are 

1. 0(/i, /2) = ^1 + h- Then the total mass of /i^ coincides with the total length of edges 
emitted in G(P) by points in V H Qx- Clearly, the so-defined is homogeneous of 
order 1 and thus Remark 1 applies. 

2. More generally, 4>{li,l2) = (h + ^2)", a > 0. Again, the total mass of /i^ is seen 
here to be the sum of power-weighted lengths of edges emitted in G{V) by points in 
V n Qx. The so-defined is homogeneous of order a. 

3. 0(/i,/2) = ^{ii+i2>e}, where 9 is some fixed threshold parameter. In this set-up, the 
total mass of /i^ is the number of edges in G{V) emitted from points in P fl Qx and 
of lengths exceeding threshold 9. This is not a homogeneous functional. 

The main tool used in our argument below is the concept of stabilization expressing in 
geometric terms the property of rapid decay of dependencies enjoyed by the functionals 
considered. The formal definition of this notion and the proof that it holds for Gilbert 
tessellations are given in Section 2 below. Next, in Section 3 the proofs of our Theorems 
1, 2 and 3 are given. 



2 Stabilization property for Gilbert tessellations 

2.1 Concept of stabilization 

Consider a generic real-valued translation-invariant geometric functional ^ defined on pairs 
{x, X) for finite point configurations A' C M^ and with x & X . For notational convenience 
we extend this definition ioi x ^ X as well, by putting ^{x,X) := C,{x,X U {x}) then. 
More generally, ^ can also depend on i.i.d. marks attached to points of X, in which case 
the marked version of X is denoted by X. 

For an input i.i.d. marked point process "P on M^, in this paper always taken to be 
homogeneous Poisson of intensity r, we say that the functional ^ stabilizes at a; G M^ on 
input V iff there exists an a.s. finite random variable R[x,V] with the property that 

^{x, V n B{x, R[x, V])) = ^{x, {V n B{x, R[x, V])) U A) (9) 

for each finite A C B{x, R[x, V]Y, with A standing for its marked version and with B{x, R) 
denoting ball of radius R centered at x. Note that here and henceforth we abuse the notation 
and refer to intersections of marked point sets with domains in the plane - these are to 
be understood as consisting of those marked points whose spatial locations fall into the 
domain considered. When (9) holds, we say that R[x, V] is a stabilization radius for V at 
X. By translation invariance we see that if ^ stabilizes at one point, it stabilizes at all points 
of M^, in which case we say that C, stabilizes on (marked) point process V. In addition, we 
say that ^ stabilizes exponentially on input V with rate C > iff there exists a constant 
M > such that 

P{R[x,V]>r}^Me~^'' (10) 

for all a; G M^ and r > 0. Stabilizing functionals are ubiquitous in geometric probability, 
we refer the reader to [1, 13, 14, 15, 16, 17, 18, 19] for further details, where prominent ex- 
amples are discussed including random geometric graphs (nearest neighbor graphs, sphere 
of influence graphs, Delaunay graphs), random sequential packing and variants thereof. 
Boolean models and functionals thereof, as well as many others. 

2.2 Finite input Gilbert tessellations 

Let A' C M^ be a finite point set in the plane. As already mentioned in the intro- 
duction, each X G A' is independently marked with a unit length random vector a^ = 
[cos(q;^), sin(a;a;)] making a uniformly distributed angle a^ G [0, tt) with the x-axis and the 
so marked configuration is denoted hj X. In order to formally define the Gilbert tessellation 
G{X) as already informally presented above, we consider an auxiliary partial tessellation 
mapping G{X) : M+ -^ J-'(M^) where J-'(]R^) is the space of closed sets in M^ and where, 
roughly speaking, G{X){t) is to be interpreted as the portion of tessellation G{X), identi- 
fied with the set of its edges, constructed by the time t in the course of the construction 
sketched above. 



Figure 1 




Finite input Gilbert tessellation. 



We proceed as follows. For each x = {x, a^) G A? at the time moment the point 
X emits in directions d^ and —dx two segments, referred to as the x~^- and a;~-branches 
respectively. Each branch keeps growing with constant rate 1 in its fixed direction until 
it meets on its way another branch already present, in which case we say it gets blocked, 
and it stops growing thereupon. The moment when this happen is called the collision 
time. For t ^ by G{X)(t) we denote the union of all branches as grown by the time t. 
Note that, with X = {xi, . . . , Xm}, the overall number of collisions admits a trivial bound 
given by the number of all intersection points of the family of straight lines {{xj + sdj s G 
M}; j = 1,2, .. . ,m} which is m{m — l)/2. Thus, eventually there are no more collisions 
and all growth unfolds linearly. It is clear from the definition that G{X){s) C G{X)(t) for 
s < t. The limit set G{X){+oo) = IJteM ^('^)(^) i^ denoted by G{X) and referred to as 
the Gilbert tessellation. Obviously, since the number of collisions is finite, the so-defined 
G{X) is a closed set arising as a finite union of (possibly infinite) linear segments. For 
X G A? by ^^{x, X) we denote the length of the upper branch x~^ emanating from x and, 
likewise, we write ^~{x,X) for the length of the corresponding lower branch. 

For future reference it is convenient to consider for each x & X the branch history 
functions x^{-), x~(-) defined by requiring that x'^^^^t) be the growth tip of the respective 
branch x'^^" at the time t G M+. Thus, prior to any collision in the system, we have just 



X 



+/- 



(t) = x + / — dxt, that is to say all branches grow linearly with their respective speeds 
+/ — dx. Next, when some y^^~ , y E X gets blocked by some other x~^^~ , x E X at time 
t, i.e. y^^^{t) = x^^~(s) for some s ^ t, the blocked branch stops growing and its growth 
tip remains immobile ever since. Eventually, after all collisions have occured, the branches 
not yet blocked continue growing linearly to oo. 
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2.3 Stabilization for Gilbert tessellations 

We are now in a position to argue that the functionals ^+ and ^^ arising in Gilbert 
tessellation are exponentially stabilizing on Poisson input V = Vr with i.i.d. marking 
according to the usual rules. The following is the main theorem of this subsection. 

Theorem 4 The functionals ^+ and ^~ stabilize exponentially on input V. 

Before proceeding to the proof of Theorem 4 we formulate some auxiliary lemmas. 

Lemma 1 Let X he a finite point set in M? and X the marked version thereof, according 
to the usual rules. Further, let y ^ X. Then for any t ^ we have 

G{X){t)AGiXU{y})it)cBiy,t) 

with A standing for the symmetric difference. 

Proof For a point set 3^ C M^ and a; G 3^ we will use the notation {x, y)^ and {x, y)^ 
to denote, respectively, the upper and lower branch outgrowing from x in G{y). Also, we 
use the standard extension of this notation for branch-history functions. Note first that, 
by the construction of G{y) and by the triangle inequality 

{x,yns') e B{y,s')^\/,^A^,yns) e B{y,s), s'^0,ee {-1,+1}. (11) 

This is a formal version of the obvious statement that, regardless of the collisions, each 
branch grows with speed at most one throughout its entire history. 

Next, write X' = X U {y} and A(t) = G{X){t)AG{X'){t) for t ^ 0. Further, let 
ti < t2 < ^3 < • • • < ^n be the joint collection of collision times for configurations X and 
X'. 

Choose arbitrary p G A(t). Then there exist unique 3^ = y{p) G {X, X'} and a; G 3^ as 
well as £ G {+, — } with the property that p = {x,yy{u) for some u ^ t. We also write 
y for the second element of {X, X'}, i.e. {3^, 3^'} = {X, X'}. With this notation, there is 
a unique i = i{p) with tj marking the collision time in y' where the branch {x, y'Y gets 
blocked in G{y'), clearly u > ti then and for s < tj we have (x, yY{s) ^ A(t). 

We should show that p G B{y,t). We proceed inductively with respect to i. For i = 
we have x = y and y = X' . Since {%j,X'Y{fS) = y G B{y,0), the observation (11) implies 
that p = {y,Xy{u) G B{y,u) C B{y,t). Further, consider the case i > and assume 
with no loss of generality that y{p) = X, the argument in the converse case being fully 
symmetric. The fact that p G G{X)(t)AG{X')(t) and that p = {x,XY{u) implies the 
existence of a point z & X' such that a branch emitted from z does block x^ in G{X') 
(by definition necessarily at the time ti) but does not block it in G{X). In particular, 
we see that {z,X')^{s) = {x,XY(ti) and {z,X')^{s') G A(s') for some 6,s,s' such that 
6 G {+, — } and s' < s ^ tj. By the inductive hypothesis we get {z,X'Y{s') G B{y,s'). 
Using again observation (11) we conclude thus that (x, A')^(tj) = {z,X'Y{s) G B{y,s) and 
hence p = {x,XY{u) G B{y,u) C B{y,t). This shows that p G B{y,t) as required. Since 
p was chosen arbitrary, this completes the proof of the lemma. □ 

Our second auxiliary lemma is 



Lemma 2 For arbitrary finite point configuration X CM."^ and x E X we have 

t{x, X) = t{x, X n B{x, 2e+(x, X))) 
C{x, X) = r(^, X n 5(x, 2i-{x, X))). (12) 

Proof We only show the first equahty in (12), the proof of the second one being fully 
analogous. Define A{X,x) = X\ B{x^2^^{x,X)). Clearly, A{X,x) is finite and we will 
proceed by induction in its cardinality. 

If |A(A?,a;)| = 0, our claim is trivial. Assume now that |A(A?,a;)| = n for some n ^ 1 
and let y = {y,ay) G A{X,x). Put t = ^~^{x,X) and X' = X\{y}. Applying Lemma 1 
we see that G{X){t)AG{X'){t) C B{y,t)._We claim that ^+{x,X) = ^+{x,X'). Assume 
by contradiction that ^^{x,X) ^ ^^{x,X'). Then for arbitrarily small e > we have 
{G{X){t)AG{X'){t))nB{x,t+e) ^ $. On the other hand, since ||x-i/|| > 2t asy ^ B{x,2t), 
for ^0 > small enough we get B{x, t + cq) fl B{y, t) = 0. Thus, we are led to 

^ {G{X){t)AG{X'){t)) n B{x,t + eo) C B{y,t) f] B{x,t + eo) = ^ 

which is a contradiction. Consequently, we conclude that t = ^^{x,X) = ^^{x,X') as 
required. Since|A(A",a;)| = n — 1, the inductive hypothesis yields ^'^{x,X') = ^^{x,X' fl 
B{x, 2^+(x, X')) = ^+(x, X' n B{x, 2t)). Moreover, X' n B{x, 2t) = X n B{x, 2t). Putting 
these together we obtain 

e{x, X) = e+(x, X') = e+(x, X' n B{x, 2t)) = e+(x, X n B{x, 2t)), 

which completes the proof. □ 

In full analogy to Lemma 2 we obtain 

Lemma 3 For a finite point configuration X dM? and x E X we have 

^+{x,X)=^+{x,XUAi) and C{x,X) = C{x,X U A2) 

for arbitrary Ai C B{x,2^+{x,X)Y, A2 C B{x,2^-{x,X)Y. 

Combining Lemmas 2 and 3 we conclude 

Corollary 1 Assume that finite marked configurations X and y coincide on B{x, 2^+(a;, X)). 
Then 

t{x, X n B{x, 2e+(a;, X))) = t{x, X) = t{x, y). 

Analogous relations hold for ^~ . 

We are now ready to proceed with the proof of Theorem 4. 



Proof of Theorem 4 We are going to show that the functional C,^ stabilizes expo- 
nentially on input process P. The corresponding statement for C,^ follows in full analogy. 
Consider auxiliary random variables C,t, Q > given by 

^t = e{x,VnB{x,g)) 

which is clearly well defined in view of the a.s. finiteness of "P fl B{x, g). We claim that 
there exist constants M, C > such that for g ^ t ^ 



P(e > t) ^ Me 
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(13) 
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Figure 2 

Indeed, let ^ ^ 0. Consider the branch a;"*" := {x^V fl B{x^ g))^ and planar regions Bi and 
Dj, i ^ 1 along the branch as represented in figure 2. Say that the event £i occurs iff 

• the region Bi contains exactly one point y oi V and the angular mark ay lies within 

(a^ + 7r/2 - e, a^c + 7r/2 + e), 

• and there are no further points of V falling into Di, 

where e is chosen small enough so as to ensure that with probability one on Si the branch 
x~^ does not extend past B^, either getting blocked in Bi or in an earlier stage of its growth, 
for instance e = 0.01 will do. Let p stand for the common positive value of P{Si), i ^ 0. 
By standard properties of Poisson point process the events Si are collectively independent. 
We conclude that, for N 9 n ^ g/3 



p(e+>3n)^p fi^n =(i-pr 



\i=l 



which decays exponentially whence the desired relation (13) follows. 

Our next step is to define a random variable i?^ = -R^[x, P,/i] and to show it is a 
stabilization radius for C,~^ at x for input process V. We shall also establish exponential 
decay of tails of R^. For ^ > we put i?^ = 2^+. Further, we set g = inf {m G N | R^ ^ m}. 
Since P{f]m(zn{^m ^ ""^D ^ infmeNP(-Rm ^ ''^) which is by (13), we readily conclude 
that so defined g is a.s. finite. Take 

i?+ := Rp (14) 

Then, using that by definition i?"*" ^ g, for any finite A C B{x, R'^Y we get a.s. by Lemma 
3 and Corollary 1 

^+{x, {V n B{x, i?+)) U A) = ^+(x, P n B{x, g) n 5(a;, 2^+(x, P n 5(a;, ^))) U A) = 

= ^+(x, P n B{x, g) n 5(a;, 2^+(x, V n fi(a;, ^)))) = ^+(x, (P n 5(a;, i?+))). 

Thus, R^ is a stabilization radius for ^+ on P as required. Further, taking into account 
that R'l = i?^ for all k ^ g by Corollary 1, we have for m G N 

P(i?+ ^ m) = P( lim Rl^m)= lim P(/2+ ^ m) = 

fc— ^-oo fe— >oo 

= lim P(e^ ^ m/2) ^ Me-^'"/^ (^5) 

whence the desired exponential stabilization follows. □ 

Using the just proved stabilization property of ^'^ and ^'^ we can now define 

^+(x, V) = ^+{x, V n B{x, R+)) = lim ^+{x, V n B{x, g)) = R+ /2 (16) 

(3— >00 

and likewise for ^~ . Clearly, the knowledge of these infinite volume functionals allows us 
to define the whole-plane Gilbert tessellation G{V). 

3 Completing proofs 

Theorems 1,2 and 3 are now an easy consequence of the exponential stabilization Theorem 
4. Indeed, observe first that, by (1), (16) and (15) the geometric functional 

ax,X):=^{e{x,X),r{x,X)) 

satisfies the p-th bounded moment condition [19, (4.6)] for all p > 0. Hence, Theorem 1 
follows by Theorem 4.1 in [19]. Further, Theorem 2 follows by Theorem 4.2 in [19]. Finally, 
Theorem 3 follows by Theorem 4.3 in [19] and Theorem 2.2 and Lemma 4.4 in [13]. 
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